Interface energies in Ising spin glasses 
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The replica method has been used to calculate the interface free energy associated with the change 
from periodic to anti-periodic boundary conditions in finite-dimensional spin glasses. At mean-field 
level the interface free energy vanishes but after allowing for fluctuation effects, a non-zero interface 
free energy is obtained which is significantly different from numerical expectations. 
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A central concept in the droplet picture of spin glasses 
is the interface free energy QjQjSSIII 5F, and the asso- 
ciated stiffness exponent 9 defined by 6F ~ where I is 
the length scale of the excitation. If 6* > the spin glass 
state is stable at finite temperature, whereas if < 
at T = large scale excitations cost very little energy so 
the spin glass state will be unstable at finite temperature. 
The value of at T = has been estimated numerically, 
in many calculations, for short range spin glass models 
from the effects of changes in boundary conditions, see 
e.g. Refs. I^SBI^- It is therefore surprising that no 
attempt has so far been made to determine the interface 
free energy from boundary condition changes using the 
alternative "Replica Symmetry Breaking" (RSB) [lOllll| 
scenario for the spin glass state. In this paper we calcu- 
late 9 analytically in high dimensions using the replica 
method and show that it conflicts with expectations from 
the droplet picture, and numerical work on systems in 
lower dimensions. 

We define the interface free energy in the standard way 
as the root mean square change in the free energy of a 
spin glass when the boundary conditions along one direc- 
tion (the z direction) are changed from periodic to anti- 



the logarithm, we have 



periodic, i.e. 5F — y AFp (here and in the following, 
the overbar means averaging over bond configurations) 
where AFp^ap = Fp — Fap, and Fp and Fap are the 
free energies with periodic and anti-periodic boundary 
conditions respectively. Anti-periodic boundary condi- 
tions can be realized by reversing the sign of the bonds 
crossing a plane whose normal is parallel to the given 
direction. It follows that AFp ap = 0. We note that ear- 
lier attempts to calculate a defect energy 0, 0| did not 
employ a definition of it which is relevant to the droplet 
picture or numerical studies. 

It is convenient to replicate the system with periodic 
boundary conditions n times and the system with anti- 
periodic boundary conditions m times, and keep n dis- 
tinct from m for the time being. Expanding the repli- 
cated partition function in powers of m and n, and taking 



{n + m) I3F 
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where AF^ = Fp — Fp = F'^p — Fap is the (mean 
square) sample-to-sample fluctuation of the free energy, 
the same for both sets of boundary conditions P or AP, 
and F = Fp = Fap- Hence, to find the variance of the 
interface free energy, AFp^p, we expand out InZ'pZ'^p 
to second order in the numbers of replicas, n and m, sepa- 
rate out the pieces involving the total number of replicas 
n + m, and take the remaining piece, which is propor- 
tional to nm. 

Using the standard replica field theory we can 
write ZpZ'^p — J Vq exp(-/3Hrop) where "Hrop is the 
replica free energy, expressed in terms of the spin glass 
order parameter field, qajsix). It is given by 
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where gap is a symmetric matrix with qaa = 0, we have 
omitted some unimportant terms of order g^, and set 
T = 1 — T/Tc- The fourth order term included is the 
one responsible for replica symmetry breaking. The co- 
efficients w and y are arbitrary positive parameters. The 
replica indices go a, 7 = 1, 2, • • • ,n,n + 1, ■ ■ ■ ,n + m. 
The order parameter q divides naturally into blocks of 
size n and m. From now on, Greek indices will label the 
first block, Roman ones the second block, so, for example, 
qai, means a € [1, and i S [n + l,n + m], and refers to 
the respective entry in the off-diagonal, or mixed, sector. 

We shall assume that there is only spatial variation in 
the z direction, which we shall take to be of length L. All 
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directions perpendicular to the z direction are of length 
M. The volume of the system is V = W^'^L. Along 
the z-direction, we impose the boundary condition that 
the solution is periodic in the Greek and Roman sectors, 
and is antiperiodic in the mixed sectors reflecting the 
sign reversal of the plane of bonds in the one sector with 
respect to the other: 



<l^jiz) = qij{z + L) 

qaiiz) = -qaiiz + L). 



(3) 



At mean-field level, there is the following stable solu- 
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tion for InZ'^Z'^p 
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is independent of the spatial coordinates and Q*^") is a 
Parisi symmetry broken saddle point solution of size sx s, 
with the necessary modification for finite positive s as 
derived in J^, i.e. 



= lim 

p — *oo 



/ p(s/p) 
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where P'^'^/'p') is a 'standard' Parisi matrix. The limit 
p ^ 00 in Eq. © should be interpreted in the same sense 
as for a standard replica symmetry breaking procedure, 
i.e. as taking p to infinity when it is convenient during a 
calculation. 

It is natural that the diagonal blocks are the same as 
the regular Parisi ansatz because ordering in the sys- 
tem with periodic boundary conditions, say, should not 
be affected by there being another completely indepen- 
dent copy with different boundary conditions. Choosing 
the mixed sector to vanish seems to be consistent with 
the standard interpretation of RSB in short-range 
systems, namely that changing the boundary conditions 
changes the system everywhere. More precisely the sur- 
face of the domain wall separating the regions which flip 
from the regions which don't flip is space filling. In this 
situation, one can reasonably expect zero overlap be- 
tween configurations with different boundary conditions. 

This solution is identical to the solution one obtains 
using the correct way of breaking the symmetry, as pre- 
sented in , for a.n-\- m-times replicated system {n + m 
being finite) without boundary condition chan ges . We 
can therefore immediately use the result from |l5l that 
on mean- field level, there is no term of order (n + m)^, let 



alone of order nm, and thus the interface energy vanishes 
to this order. 

We now turn to the loop expansion about the sad- 
dle point, which is expected to be valid for dimension d 
greater than 6. The first correction is due to Gaussian 
fluctuations around the saddle point solution. They are 
given by 
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where 



/(A:^) = ^d^ln(fc2 + A, 
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/c is a d-dimensional wave vector and A^, are the eigen- 
values of the Hessian, evaluated at the saddle point so- 
lution, and their degeneracies. The eigenvalues and 
degeneracies df_i are the same as for a system of size n + m 
without boundary condition changes (because the saddle 
point solution is the same), only the nature of the k- 
vectors changes for the terms involving eigenvalues whose 
corresponding eigenvectors / are nonzero exclusively in 
the mixed sector (i.e. fap = fij = 0): the wave vectors 
have to respect the imposed boundary conditions, which 
implies k = {2nnr/M, . . . ,2nd-in/M, {2nd + 1)7:/ L) 
(with ni G Z) in the mixed sector as opposed to fc = 
{2nin/M, 2n<i_i7r/M, 2?i<j7r/L) in the Greek or Ro- 
man sectors. 

It was shown in .15.] for a system without boundary 
condition changes that it is initially easier to compute 
dl/d{k^) than / itself , and that it is given in terms of 
the diagonal propagators Gap,ap (or Gff in the hmit of 
infinitely many replica symmetry breaking steps |l4| ) as 
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where we have dropped the subscript 11 from the prop- 
agators as it is irrelevant here and replaced it by n since 
the propagators depend on it. 

Therefore the contribution to dl/d{k^) from those 
eigenvectors that are nonzero in the Greek or Roman sec- 
tors (the periodic sectors, hence the subscript P below) 
is 



dip 

a(fc2 



dx G^ 



n 



m 



dxGl'' 



(10) 
(11) 



The last line follows from the modified symmetry break- 
ing procedure (Eq. ©), as was shown in The origin 
of the term linear in n -|- m in Eq. (|10|l is the eigenvectors 
that are nonzero in a Parisi block p("/p) or p('"/p) on the 
diagonal while the origin of the n^ + m^-term is the 
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eigenvectors that are nonzero in the ofF-diagonal blocks. 
This observation faciUtates calculating the contribution 
from the mixed sector as there are only eigenvectors of 
the latter type present, i.e. there is no term of linear or- 
der. Therefore dlAp/d{k^) is given by 



dlAP 



(12) 



where the prefactor nm reflects the number of eigenvec- 
tors in the mixed sector. 

The integral / d{k'^)G^ and the constant of integra- 
tion have been worked out in |l5j , resulting in 



J{k^) := / rf(fc2)GS" = Hk^ + ^) 

: tan 



are in fact dominating over the terms we kept. Therefore 
we need to treat the fc = mode properly before pro- 
ceeding. The way to do this is to go to the equation of 
state for qaj3 and include the A; = mode exactly, while 
treating the other modes perturbatively as before. The 
complete equation of state is given by Eq. (15) from [l^. 
and restricted to the k — Q mode it reads 



WXi 



(13) 

where Xi is the breakpoint of the Parisi g-function. We 
can now assemble in / the terms of quadratic order. 



/ (n + m)C + 
^{n + m)C + 



Jp{k^) + nmJAp{k^) (14) 



[n + m)^ 



Jp(fc2) + nm{JAp{k^) - Jp{k^)) 

(15) 



The constant C is of no interest to us. The subscripts 
P and AP on J mean that J must be taken as when 
the argument is not of the required type, i.e. periodic or 
antiperiodic. 

We can now identify the term that gives rise to the 
interface energy. Comparison with Eq. shows 
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where the subscripts on the sums indicate the nature of 
the allowed fc- vectors, as made explicit in the second part 
of the equation where the z component of the /c-vector 
has been split off, leaving the d — 1-dimensional wave vec- 
tor I. The sum over the z component has been extended 
to ±oo, introducing only exponentially small errors for 
large L. 

We note a potential pitfall in this result: the contribu- 
tion to Eq. Hlt)|) from the fc = term (in X]p) diverges. 
Usually, this problem is removed by converting the sums 
to integrals converging in high enough dimensions, and 
arguing that the divergence is, in reality, only a subdom- 
inant contribution. However, since B < {d~l)/2 4], the 
interface energy is subdominant itself, so it is not clear 
whether the subdominant terms from the fc = mode 



I ^ E '^a-t,l3-t{k = 0) 4- 2yqai3Gaf3,af3{k = 0) 



V 



(17) 

This equation is highly nontrivial since G in this expres- 
sion is the full propagator. We do not propose to solve 
this formidable self-consistency equation, but we note 
that the presence of the right hand side shifts Qap by 
an amount eafj from the mean-field value, which in turn 
shifts the eigenvalues of the Hessian. The left hand side 
is given by T^-yS^'p^^si'^ = 0)e^5 = 0{e) (recalling that 
G~^{k = 0) is equal to the Hessian), the right hand side 
is of order l/FAmin, where Amin is the smallest eigen- 
value of the Hessian. If Amin — 0{e), which is the natu- 
ral expectation, it follows that Amin ~ V"^^'^. Therefore 
G(fc — 0) has changed from being infinite to being of or- 
der xhis argument is not rigorous, however, there- 
fore we prefer to denote the exponent more generally by 
2fj,. The upshot of this treatment is that we can exclude 
the divergent fc = terms from the sums over wave vec- 
tors (as they have been dealt with non-perturbatively), 
provided a term of order V'^^, where /x may be 1/4, is 
introduced in the and terms in Eq. (|14|l . This ad- 
ditional term is identical to the free energy fluctuations 
in the Sherrington-Kirkpatrick model which has only the 
fc = mode, and will be denoted by Afg^V'^'^. This ob- 
servation allows us to obtain estimates of /i from existing 
numerical work 0,0,0], which supports fi — 1/4. 

Since we are expecting that the changes to the eigenval- 
ues are of order V~^^^, while the changes due to the dif- 
ferent boundary conditions are of order 1/ L'^, our treat- 
ment of the non-zero fc modes will be satisfactory in the 
range of dimensions where the loop expansion applies, 
i.e. d > 6. 

Upon completing the square as in Eq. H15I) the con- 
tribution A/Ij^U^m appears in the nm term, so from 
Eq. Jini) we get 



P,AP 

oc 
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The sums over r in Eq. lfTS|l can be calculated exactly, 
in principle, and the sum over / can be converted to an 
integral with a lower cutoff and carried out, but the re- 
sult is too long to show here. The important feature 
of it is that the leading behaviour as a function of L is 
determined by the divergent part of J as fc^ ^ 0. The 
other parts of J only give exponentially small corrections. 
Since J(fc^) « —Tiw/Ayk'^ for small fc^, it is sufficient to 
work out the term 



Ay 



E 



1 



1 



wnL 



(2r+l)2 



Together with J27Li^ 



L2 



(2r-l)^ 



/2 + 



L2 



iyl sinh IL 
(19) 
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7r^/12 this gives 



d-l 



dll 



d-2 



wnL 



Ayl sinh/i 



24^ 



(20) 



iV^(L/M)-^i2 + A/|KF^^ (21) 



where 



f{L/M) = 



Ay{2TT)d~^ 



dx : 



sinhx 



(22) 



is an exponentially decreasing scaling function and Sd is 
the surface of a d-dimensional unit sphere. 

Only the first term in Eq. I|21l) has a form compat- 
ible with aspect ratio scaling 20], according to which 

^JAF^^^p = L^fiL/M). On the face of it, this would 
give rise to 9 — 1 for all dimensions. The other two terms, 
however, do not have aspect ratio scaling form. In par- 
ticular, the term A/gj^F^^, which is dominant in d > 6 
if /Lt > 1/6, depends only on volume but not on shape. 

Our calculation is exact in high dimensions within the 
replica symmetry breaking scenario for spin glasses. It is 
quite unusual and contradicts all expectations one might 
have about the interface energy based on experience from 
other systems and numerical data. It is significantly dif- 
ferent from that found in, for instance, ferromagnets. 
There, the defect energy comes from the gradient term 
in the analogue of Eq. ||2Jl. Here, on the other hand, the 
mean-field solution is independent of z so there is no con- 
tribution from the gradient term. A difference between 
the interface energy in spin glasses and ferromagnets is, 
however, that in ferromagnets there is a 'real' domain 
wall, whereas in spin glasses, the interface can only be de- 
fined by comparing one system to a reference system with 
the opposite set of boundary conditions. Thus strictly 
speaking, the interface in spin glasses is not a physical 
system itself, which may account for the absence of an 
interface energy on the mean field level. 

The failure of aspect ratio scaling is a strong prediction 
which contradicts numerical evidence for d = 2 12(1 |2l| . 



The replica symmetry breaking scenario predicts space- 
filling domain walls [la, |22| , therefore the dependence of 
the interface energy on volume but not on shape (to lead- 
ing order) appears natural since the interface explores 
even the remote corners of the sample and would be 
likely to exist in some form down to three dimensions, 
even though the loop expansion used in this paper will 
need modification below six dimensions. This suggests a 
simple test of replica symmetry breaking ideas. If they 
are valid in three dimensions, then aspect ratio scaling 
will fail. To date, there is (weak) evidence that aspect 
ratio scaling works in three dimensions |20| . 
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